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ABSTRACT 

Core collapse of dense massive star clusters is unavoidable and this leads to the 
formation of massive objects, with a mass up to 1000 M & and even larger. When 
these objects become stars, stellar wind mass loss determines their evolution and 
final fate, and decides upon whether they form black holes (with normal mass 
or with intermediate mass) or explode as a pair instability supernova. In the 
present paper, we discuss the evolution of very massive stars and we present a 
convenient evolution recipe that can be implemented in a gravitational N-body 
code to study the dynamics of dense massive clusters. 

Subject headings: stellar evolution - stellar wind - very massive stars 



1. Introduction 

The inner 100 pc of the Galactic center contains several yo ung dense star cluster s 



(IFiger et al.l Il999al ). some of them with reliab le mass estimate s ( IBorissova et al.l 120051 ). 



Of particular inter est are th e Arches cluster (jFiger et al.l 120021 ). the Quintuplet cluster 



(IFiger et al.lll999bh . IRS 13E flMaillard et al. 1120041 ) and IRS 16SW flLu et al. 1120051 ). Grav- 



itational N-body simulations reveal that soon after birth such clusters may experience 



stars participate in a "collision runaway" or "collision merger" 


(Qu 


inlan & Shapiro 


1990 


Porteeies Zwart et al. 1999; Figer et al. 


2002: 


Gurkan et al. 


2004 




Gurkan & Rasio 


2005; 


Freitag et al. 


2006 


). In a recent paper, '. 


3 ortegies Zwart et al. 


(• 


zooe 


) estimated the typical 



mass of these objects. They concluded that clusters in the inner 10 pc (respectively between 
10 pc and 100 pc) of the Galactic center form collision mergers with an average mass ~ 1000 
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M & (respectively ~ 500 M ). The dynamical evolution of a cluster where core collapse hap- 
pens obviously depends on whether or not the very massive merger becomes a very massive 
star and, when a very massive star is formed, on the evolution of this very massive star. 

The present paper deals with the evolution of very massive stars, products of runaway 
merging. The computation method is outlined in section 2 where we provide an easy evolu- 
tion calculation recipe. It is obvious that this evolution is critically affected by stellar wind 
mass loss. The formalism that we use is discussed in section 3. The evolution of stars with 
a post-merger mass between 300 M and 1000 M & is illustrated in section 4. 



2. Simulating the evolution of very massive stars 

In the subphotospheric layers of very massive stars, where the opacity becomes larger 
than the electron scattering value, the radiation force almost balances gravity, causing a 
core/extended halo stellar structure. This hampers the convergence of stellar evolutionary 
computations. However, since the mass of these layers is very small, they hardly affect the 
overall internal structure and treating these layers using Thomson scattering opacity only, 
still provides an accurate description of very massive stellar evolution while at the same time 
avoiding any numerical difficulties. 

Our calculations of t he ey olution of very massive stars are based on the results of 



Nadyozhin fc Razinkova I (120051 ) (NR) who constructed interior models for massive objects 
using the similarity theory of stellar structure (treated as a boundary- value problem). Their 
models correspond to chemically homogeneous stars, having Thomson scattering as the only 
opacity source throughout. The obtained model sequences depend on one parameter only; 
/i 2 M (/i being the mean molecular mass of the gas and M the total mass of the star). 
During most of their evolution very massive stars produce a convective core that almost 
covers the entire star, meaning that their evolution can be simulated accurately with a 
homogeneous model. The fact that very massive stars are expected to lose a significant 
amount of mass by stellar wind (section 3) strengthens the conclusion that very massive stars 
evolve in a quasi- homogeneous way. Furthermore, in case of steep dependencies of the energy 
generation rate on temperature (as is the case for the CNO cycle and the 3a-reaction), most 
of the energy production is localised near the very centre of the star. Therefore, to a very 
good approximation, the luminosity is constant throughout the star and the dimensionless 
luminosity equation is decoupled from the rest of the stellar structure equations. This means 
that the model sequence of NR can be used to describe the core hydrogen burning (CHB) 
as well as the core helium burning (CHeB) stage of very massive stars. 
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To simulate the evolution of a star, we proceed as follows. NR provide best fit relations 
(see their eqs. 30, and 34 in combination with 36) as a function of fi 2 M for (a.o) the 
luminosity and convective core mass of their computed sequence (which covers the range 
< fi 2 M < 4000 Mq). For a given stellar luminosity (and assuming that central nuclear 
burning is the only energy source in the star), conservation of energy allows one to derive 
the amount of nuclear fuel that is burned per unit of time. Then, from knowledge of the size 
of the convective core, a differential equation for the variation of the central abundance of 
the fuel is obtained. For the CHB case, we have: 



M cc (/i,M) — = (1) 

at €h 

where e# is the energy produced from burning one mass unit of hydrogen. In this equation, 
both M cc and L vary due to changes not only of (/,, but also of M through stellar wind mass 
loss: 



dM 
~~db 



M(/i, M, L). 



(2) 



Assuming the mass loss rate M can be derived from knowledge of /i, M and/or L, 
the solution of this coupled set of two differential equations provides the evolution of the 
quantities X, fi, M, M cc , L and M as a function of time, as well as the duration of the 
Hydrogen burning stage (which ends, of course, when X = 0). 

For the CHeB case, one needs to account for the fact that C and O are produced in 
a non-constant ratio, which affects the energy product ion per unit ma ss of burned helium. 
The equivalent of eq. 1 for Helium burning is (see also lLanger I (ll989al ) ) : 



/ »^f B YdY B c dC B dO\ T , 

In eq. 3, the B symbols represent the binding energies of the corresponding nuclei 
whereas the A symbols are their atomic weights. We impose the additional contraint Y + 
C + O = 1 and thus dY + dC + dO = 0, for simplicity. 

Langer (1989b) computed models for massive homogeneous CHeB stars and found an 
abundance evolution of C and O relative to He that was closely followed by all models 
(15 Mq < M < 100 Mq), independent of initial mass. This resulted in a fit between C and 
Y (his eq. 1, see also his Fig. 1). According to Langer, this fit is very accurate for Y 
values larger than about 0.5, when the O abundance is low and the 12 C(a,7) 16 reaction 
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is relatively unimportant. At lower values of Y, this reaction and the fact that some 16 O is 
converted into 20 Ne towards the end of CHeB of very massive stars, produce an estimated 
uncertainty in the fit of about 5-10%. Using this C(Y) relation and eliminating the C and 
O time derivatives from eq. 3, one finally obtains: 

Here C"(F) denotes the derivative of the C(Y) fit of Langer, with respect to Y. As was 
the case for H burning, the combination of eq. 4 with a mass loss rate formalism of the form 
of eq. 2 enables one to compute the evolution of the star upto He depletion in the core. 

In section 4 we will compare evolutionary results of massive stars which are calculated 
with the similarity theory with results calculated with detailed stellar evolutionary codes in 
order to evaluate our computational method. 



3. The stellar wind mass loss formalism of very massive stars 

Since direct observations of very massive stars in general, their stellar wind mass loss 
rates in particular are lacking, we are forced to estimate the effect of stellar wind mass 
loss on very massive star evolution either by extrapolating empirical for malisms hold i ng for 



massive stars or we can use theoretical models when they are available. iKudritzki I (120021 ) 
studied line-drive winds of very massive stars and calculated mass loss rates as function of 
metallicity Z of very massive O-type stars with a luminosity Log L/L between 6.3 and 7.03 
and for T e jy- values between 40000K and 60000K and he presented M-interpolation formulae 
for three different T e //-values. Notice that the mass loss rates are smallest for the highest 
T e ff- This means that if, by using the M-interpolation formula corresponding to the highest 
T e ff during the whole CHB phase, we predict that stellar wind mass loss of very massive 
stars is large, the real mass loss may be even larger. 

Due to stellar wind mass loss during CHB, the post-CHB (CHeB) remnants of the very 
massive stars are hydrogen deficient and may be considered as very massive Wolf-Rayet (WR) 
stars, e.g. the evolution of very massive stars during CHeB has to be calc ulated accounting 



for W R-like mass loss rates. Theoretical formalisms have been presented by lNugis fc Lamers 



( 120021 ) but they contain stellar wind parameters which can only be derived by linking a full 
hydrodynamical wind model to a stellar evolutionary model (e.g. the temparature and radius 
in the wind at the sonic point, the wind velocity at infinity etc). Therefore we prefer to rely 
on empirical formalisms. Based on indirect arguments involving population synthesis of WR 
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stars in the Solar neighbourhood (the WN/WC number ratio), on the masses of black holes 
in binaries and on direct mass loss rate determinations of WR stars including the effects of 
clumping (before 1998 the effects of clump i ng on empirical mass loss rates was investigated fo r 



only a few WR star), Vanbeveren et al. ( 19981 ) (see also Van Bever &: Vanbeveren ( 20031 )) 



proposed the following relation 



Log(M) = Log(L) - 10 + 0.5Log(Z/Z Q ). 



(5) 



where Z stands for the initial metallicity which is proportional to the Fe abundance of the 
WR star. 



Nugis &: Lamers I (120001 ) used dumping-corrected mass loss loss rates of a large sample 



of Galactic WR stars and proposed the following M formula as function of luminosity and 
helium abundance Y 



Log(M) = -11 + 1.29Log(L) + 1.7 Log(Y) + 0.5Log(l — X — Y). 



(6) 



In the two formulae given above M is in M /yr and L in L & . 



Remarks. iKudritzki I (120021 ) calculated the mass loss rates for stars with a luminosity 
Log(L/ Lq) up to 7. Initially (on the zero age main sequence) our 1000 M star has 
Log(L/ L Q ) = 7.5 and we extrapolated the M-interpolation formulae. We obviously as- 
sured that the mass loss rates in the extrapolation zone are smaller than or equal to the 
maximum mass loss rates for line driving as discussed by lOwocki et al. I (120041 ) . 



The WR stars where both empirical WR mass loss rate formulae given above hold, have 
luminosities in the range 5.0 < LogL < 6.0. The theoretically predicted very massive WR 
stars (section 4) have luminosities up to Log L = 7. The very massive WR mass loss rates 
that we use here are therefore extrapolated values implying quite some uncertainty. In the 
next section we will discuss the consequences of this uncertainty. Similarly as for the CHB 
mass loss rates it is obvious that also here we check that the extrapolated values remain 
smaller than or equal to the maximum rates. 



Results 



In order to illustrate to what extent our very massive star evolutionary scheme ap- 
proaches detmkd_evcJutiOTwy computations, Table 1 compares the 120 M & evolutionary 
result of ISchaller et al. I (119921 ) with our prediction where we obviously used the same mass 
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loss prescription during CHB and during CHeB as in Schaller et al. As can be noticed, the 
correspondence is very good. The basic assumption of our method is the quasi-homogeneous 
evolution of very massive stars. The more massive a star, the larger is the convective core 
and the larger is the stellar wind mass loss rate. This means that the more massive a star 
the closer its evolution will be to the quasi-homogeneous one. Since our method gives very 
satisfactory results already for the 120 M Q star we are inclined to conclude that it will closely 
describe the e v olutio n of very massive stars. The latter is strengthened by the following. 



Marigo et al. I (120031 ) calculated the evolution of a 1000 M & zero-metallicty star which is 
subject to a large stellar wind mass loss (they also use the Kudritzki formalism). In Ta- 
ble 1 we also compare the results of Marigo et al. with ours whereas Figure 1 compares 
the temporal behaviour of the mass of the convective core, the luminosity and the effective 
temperature. As can be noticed, the correspondence is excellent. 

We calculated the evolution of stars with a mass in the range 300-1000 M for three 
metallicities, Z = 0.04, Z = 0.02 and Z = 0.001, using the core hydrogen burning M inter- 
polation formulae corresponding to the highest T e ff (section 3). The collision of massive 
stars in a dense cluster happens typically 1 or 2 Myrs after their formation. The central 
hydrogen abundance of the stars at the moment of collision may be significantly smaller than 
the initial value which implies that even when the merger product is well mixed and becomes 
homogeneous, the new X may be significantly smaller than the initial value of the cluster. 
For this reason we performed evolutionary calculations of very massive stars with X = 0.68, 
0.6 and 0.5. The results are summarized in Table 2, the initial mass - final mass relationship 
is depicted in Figure 2 (in this figure we also plot the relation for stars with initial mass 
smaller than 100 M Q , taken from Van Bever and Vanbeveren, 2003), the temporal evolution 
of the stellar mass is shown in Figure 3 and the evolutionary behaviour in a mass-luminosity 
diagram is given in Figure 4. They illustrate the following conclusions: 



• Very massive stars with Z > 0.02 and with initial mass > 300 M Q lose most of their 
mass in the form of stellar winds during CHB and CHeB. The same applies for very 
massive stars with Z = 0.001 and with initial mass > 500 M & . 

• The final masses at the end of CHeB calculated with the WR mass loss rate formula 
5 and 6 are very similar. 

• All the very massive stars with the same initial chemical composition and with an 
initial mass > 300 M Q end their life as stars with very similar final mass and their 
CHB and CHeB lifetimes are very similar. This is the reason why in Table 2 we only 
give the details for the 300 M and 1000 M Q stars. 
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A very massive star with a lower initial X has a shorter CHB lifetime, but a larger 
luminosity, thus a higher stellar wind mass loss rate. This explains why the final masses 
for very massive stars with the same initial metallicity Z hardly depend on the initial 
X. 

Very massive OB-type stars with Z > 0.02 and WR stars with Z > 0.001 obey a very 
tight mass luminosity relation, i.e. 

Log(L) = 1.07(Log(M)) 2 - A.62Log(M) + 11.8 Z = 0.04 (7) 

Log(L) = l.\2{Log{M)f - 4.98Log(M) + 12.4 Z = 0.02 (8) 
for the OB stars and 

Log(L) = 1.23Log(M) +4.19 Z = 0.04 (9) 

Log(L) = l.ULog(M) + 4.37 Z = 0.02 (10) 

Log(L) = 0.88Log(M) + 4.98 Z = 0.001 (11) 
for the WR stars. All relations have a determination coefficient R 2 > 0.99. 

When Z is larger than or equal to 0.02 our computations reveal that the very massive 
stars will end their life as a < 40-50 M Q black hole. Since the Galactic bulge has such 
a large Z, intermediate mass black holes with a mass of a few 100 M & may be difficult 
to form there. 

For Z = 0.001 the final mass of very massive stars < 170 M & . Therefore, intermediate 
mass black holes (but with a mass of a few 100 M ) may form in dense metallicity 
poor clusters. 



From the results of iHeger fc Woosley I (120021 ) we conclude that when Z is between 



0.001 and 0.02, one may expect pair-instability supernova candidates among collision 
runaway mergers in clusters in the Galactic center. 



Remarks. The evolutionary computations discussed above let us conclude that it may be 
difficult to form intermediate mass black holes by means of runaway merging in dense clusters 
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in the Galactic bulge (where Z > 0.02). Of course the computations rely on the adopted 
mass loss rate formalisms. During CHB we used the theoretically calculated M-interpolation 
formula corresponding to the highest T e jj-values which means that in reality the overall 



mass t hat is lost during CHB may be larger than the values in Table 2. Note that iKudritzki 



(120021 ) calculated the mass loss rates of radiaton driven stellar winds. Additional processes 
(like rotation for example) will increase the derived rates. Furthermore, very massive stars 
may experience a luminous blue variable (LBV) phase somewhere near the end of CHB, 
much like the ma s sive s tars do. The LBV phase is characterised by eruptive (explosive, 



Smith fc Owocki I (120061 )) mass loss episodes (as observed in rj Car) and this may increase 
the total mass loss. The remarks discussed above stengthen the main thesis of the present 
paper. 

A major uncertainty is obviously the (extrapolated) empirical WR mass loss rate formal- 
ism for the very massive stars. To illustrate the importance of this uncertainty we computed 
the evolution of the very massive stars (Z > 0.02) but with WR mass loss rates which are 
a factor 2 and 4 smaller than predicted by the relations given in section 3. The results are 
given in Table 2 as well. As expected, the final CHeB masses are larger and some of them fall 
in the range where we expect pair-instability supernovae to happen, e.g. these stars do not 
form BHs at all. Since a pair-instability supernova happens roughly when the final CHeB 
mass is larger than 65-75 M Q , we conclude from our computations and the argumentation 
above that it is very improbable that very massive stars in the Galactic center produce black 
holes with a mass larger than 65-75 M Q . 



5. Summary 

In the present paper we studied the quasi- homogeneous evolution of very massive stars 
with mass up to 1000 M which could be the result of core collapse of young dense clusters. 
When the theory of radiatively driven stellar wind mass loss applies it follows that the 
evolution of very massive stars is dominated by these winds during CHB and during CHeB. 
At Solar metallicity and larger very massive stars end their live as a black hole with a mass 
less then 75 M Q . At Z = 0.001, the final mass of the very massive stars studied here may be 
a factor 2 to three larger compared to those at Z = 0.02, e.g. in low metallicity regions the 
formation of intermediate mass black holes with a mass of a few hundred M & is a possibility. 
Furthermore it is very plausible that between Z = 0.02 and Z = 0.001 at least some very 
massive stars will end their life with a pair-instability supernova. During a pair-instability 
supernova very large amounts of metals may be ejected and we like to suggest that the metal 
poor (Z < 0.02) chemical evolution of galactic bulges may be affected by cluster dynamics 
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through the formation of very massive stars and the occurence of pair-instability supernova. 
Acknowledgement 

We like to thank an unknown referee for very valuable suggestions that improved the scientific 
content of the paper. 

REFERENCES 

Borissova, J., Ivanov, V.D., Minniti, D., Geisler, D., Stephens, A.W.: 2005, A&A435, 95. 

Figer, D.F., Kim, S.S.: 2002, in ASP Conf. Ser. 263: Stellar Collisions, Mergers and their 
Consequences, pp. 287. 

Figer, D.F., McLea, I.S., Morris, M.: 1999b, ApJ514, 202. 

Figer, D.F., Najarro, F., Gilmore, D., Morris, M., Kim, S.S., Serabyn, E., McLean, I.S., 
Gilbert, A.M., Graham, J.R., Larkin, J.E., Levenson, N.A., Teplitz, H.I.: 2002, 
ApJ581, 258. 

Figer, D.L., Kim, S.S., Morris, M., Serabyn, E., Rich, R.M., McLea, I.S.: 1999a, ApJ525, 
750. 

Freitag, M., Giirkan, M.A., Rasio, FA.: 2006, MNRAS 368, 141. 
Giirkan, M.A., Freitag, M., Rasio, FA.: 2004, ApJ604, 632. 
Giirkan, M.A., Rasio, FA.: 2005, ApJ628, 236. 
Heger, A., Woosley, S.E.: 2002, ApJ567, 532. 
Ishii, M., Ueno, M., Kato M.: 1999, PASP 51, 417. 
Kudritzki, R.P: 2002, ApJ577, 389. 
Langer, N., A&A, 1989a, 220, 135. 
Langer N., A&A, 1989b, 210, 93. 

Lu, J.R., Ghez, A.M., Hornstein, S.D., Morris, M., Becklin, E.E., 2005, ApJ625, 51. 
Marigo,P, Chiosi,C, Kudritzki,R.-P2003, A&A 399, 617. 



-10- 

Maillard, J.P., Paumard, T., Stolovy, S,R., Rigaut, F.:,2004, A&A 423, 155. 
Nadyozhin, D.K., Razinkova, T.L., 2005, Astronomy Letters, Vol. 31, No. 10, pp. 695. 
Nugis, T., Lamers, H.J.G.L.M.: 2000, A&A 360, 227. 
Nugis, T., Lamers, H.J.G.L.M.: 2002, A&A 389, 162. 
Owocki, S.R, Gayley, K.G., Shaviv, N.J., 2004, apj, 616, 525. 

Portegies Zwart, S.F., Baumgardt, H., McMillan, S.L.W., Makino, J., Hut, P., Ebisuzaki, 
T.: 2006, ApJ(in press) ( ]astro-ph/05lT397| ). 

Portegies Zwart, S.F., Mkino, J., McMillan, S.L.W., Hut, P.: 1999, A&A 348, 117. 

Quinlan, G.D., Shapiro, S.L.: 1990, ApJ356, 483. 

Schaller, G., Schaerer, D., Meynet, G., Maeder, A.: 1992, A&AS 96, 269. 

Smith, N., Owocki, S.P: 2006, ApJ645, L45. 

Van Bever, J., Vanbeveren, D.: 2003, A&A 400, 63. 

Vanbeveren, D., De Loore, C, Van Rensbergen, W.: 1998, The A&A Rev. 9, 63. 



This preprint was prepared with the AAS IATgX macros v5.2. 



- 11 - 




t/(10° yr) 



7.5 1000 



800 



600 



400 





Similarity theory 




Marigo et al. 2003 







0.5 1 

t/(10 6 yr) 



1.5 



Fig. 1. — Left figure: the temporal behaviour of the luminosity (in L & ) and of T*ff during 



CHB of a zero metallicity 1000 M & star: a comparison between the results of iMarigo et al 
( 120031 ) (dashed lines) and our method (full lines). Right figure: the same as the left figure 
but for the total mass (upper lines) and the mass of the convective core (lower lines), (masses 
are in M@). 
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Fig. 2. — The initial mass (M init ) - final mass (M end ) relation (all masses are in M ) for Z 
= 0.04 (dotted), Z = 0.02 (full line) and Z = 0.001 (dashed line). 
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Fig. 3. — The temporal evolution of the mass (in M Q ) for the stars with initial mass = 1000 
M Q , 750 M Q , 500 M Q and 300 M . The dot-mark corresponds to the end of CHB = the 
beginning of the WR phase. 
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Fig. 4. — The evolution in the mass - luminosity diagram (both in solar units) of the stars 
with initial mass = 1000 M , 750 M , 500 M Q and 300 M . The dots give the location of 
the initial zero age parameters of the four stars. The lower right figure shows the mass - 
luminosity evolution during the He-burning phase for the 3 metallicities: Z = 0.001 (dashed 
line), Z = 0.02 (full line) and Z = 0.04 (dotted line). 
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Mi n n 


Z 


M e CHB 


TcHB 


M e wR 


TcHeB 


reference 


1000 


2.10~ 6 


466.76 


2.050 


415.18 


2.31 


This paper 






482 


1.861 


419 


2.15 


Marieo et al. ( 20031 


120 


0.02 


69.64 


2.753 


7.91 


4.65 


This paper 






70.28 


2.611 


7.77 


4.40 


Schaller et al. (1992) 



Table 1: A comparison of massive and very massive star evolutionary computations performed 
with a detailed evolutionary code and with the similarity method used in the present paper. 
All masses are in M , T C hb is in Myr, T C HeB is in 10 5 yr. 
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Mi n u 


z 


X 


M eC HB 


TcHB 


M eC HeB 


TcHeB 


1000 


0.04 


0.68 


140.09 


2.045 


22.31 - 


54.79 - 


87.12 




3.34 - 2.95 - 


2.83 








0.60 


141.64 


1.724 


22.54 








3.33 










0.50 


144.35 


1.353 


22.95 








3.31 








0.02 


0.68 


152.91 


2.017 


40.72 - 


78.37 - 


109.06 - 


35.45 


3.03 - 2.83 - 


2.76 - 


3.36 






0.60 


154.56 


1.701 


41.19 








3.02 










0.50 


157.35 


1.337 


41.86 








3.01 








0.001 


0.68 


220.93 


1.903 


163.13 








2.62 










0.60 


222.47 


1.610 


164.41 








2.61 










0.50 


225.15 


1.272 


166.42 








2.61 






300 


0.04 


0.68 


136.85 


2.131 


21.86 - 


53.61 - 


85.14 




3.36 - 2.96 - 


2.84 








0.60 


138.37 


1.786 


22.54 








3.33 










0.50 


140.95 


1.394 


22.95 








3.31 








0.02 


0.68 


148.76 


2.109 


39.77 - 


76.30 - 


106.09 - 


34.86 


3.04 - 2.84 - 


2.77 - 


3.38 






0.60 


150.20 


1.769 


40.05 








3.04 










0.50 


152.74 


1.382 


40.69 








3.03 








0.001 


0.68 


205.96 


2.038 


152.15 








2.64 










0.60 


207.30 


1.710 


153.09 








2.64 










0.50 


209.63 


1.337 


154.96 








2.63 







Table 2: Evolutionary properties of very massive stars as function of initial chemical compo- 
sition (Z, X). We list the mass at the end of CHB (M s chb), the CHB timescale (Tchb), the 
mass at the end of CHeB [M e cueB) and the CHeB timescale {TcHeB)- The CHeB numbers 
are always calculated using the WR mass loss rate formula 5. For Z = 0.04, Z = 0.02 and X 
= 0.68 we also list the two CHeB parameters using WR mass loss rate formula 5 divided by 2 
and using formula 5 divided by 4. For Z = 0.02 and X = 0.68 the fourth number corresponds 
to the case where the WR mass loss rate is computed with formula 6. All masses are in M , 
Tchb is in Myr, T CHeB is in 10 5 yr. 



